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f—N ■ In this paper, we are exploring some of the properties of the self- 

similar solutions of the first kind. In particular, we shall discuss the 
kinematic properties and also check the singularities of these solutions. 



^ | We discuss these properties both in co-moving and also in non co- 

i , moving (only in the radial direction) coordinates. Some interesting 

§j\\ features of these solutions turn up. 
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1 Introduction 

A set of field equations remains invariant under a scale transformation if we 
assume appropriate matter field. This implies the existence of scale invari- 
ant solutions to the field equations. Such solutions are called self-similar 
solutions. The special feature of self-similar solutions is that, by a suit- 
able coordinate transformations, the number of independent variables can 
be reduced by one and hence reduces the field equations. In other words, 
self-similarity refers to an invariance which simply allows the reduction of a 
system of partial differential equations to ordinary differential equations. 

Similarity solutions were first studied in General Relativity (GR) by 
Cahill and Taub [1]. They studied these solutions in the cosmological con- 
text and under the assumption of spherically symmetric distribution of a 
self-gravitating perfect fluid. They assumed that the solution was such that 
the dependent variables are essentially functions of a single independent vari- 
able constructed as a dimensionless combination of the independent variables 
and that the model contains no other dimensional constants. They showed 
that the existence of a similarity of the first kind in this situation could be 
invariantly formulated in terms of the existence of a Homothetic vector (HV). 

In GR, self-similarity is defined by the existence of a HV field. Such 
similarity is called the first kind (or homothety). There exists a natural 
generalization of homothety called Kinematic self-similarity, which is defined 
by the existence of a kinematic self-similar (KSS) vector field. The basic 
condition characterizing a manifold vector field £ as a self-similar generator 
[2] is given by 

£ ( A = \A, (1) 

where A is constant and A is independent physical field. This field can be 
scalar (e.g. //), vector (e.g. u a ), or tensor (e.g. g ab ). In GR, the gravitational 
field is represented by the metric tensor g ab , and an appropriate definition of 
geometrical self-similarity is necessary. 

A kinematic self-similarity satisfies the condition 

£^u a = au a , (2) 

with 

£ ( h ab = 25h ab , (3) 

where a and 5 are constants and h ab = g ab — u a u b is the projection tensor. 



KSS perfect fluid solutions have been explored by several authors. Carter 
and Henriksen [3] defined the other kinds of self-similarity namely second, ze- 
roth and infinite kind. In the context of kinematic self-similarity, homothety 
is considered as the first kind. 

The only barotropic equation of state which is compatible with self- 
similarity of first kind is p = kfi. Carr [4] has classified the self-similar 
perfect fluid solutions of first kind with this equation of state for the dust 
case (k = 0) and the case < k < 1 has been studied by Carr and Coley 
[5]. Coley [6] has shown that the Friedmann Robertson Walker solution is 
the only spherically symmetric homothetic perfect fluid solution in the par- 
allel case. Mcintosh [7] has discussed that a stiff fluid (k = 1) is the only 
compatible perfect fluid with the homothety in the orthogonal case. 

Benoit and Coley [8] have studied spherically symmetric spacetimes which 
admit a KSS vector of the second and zeroth kind. Sintes et al. [9] have 
considered spacetimes which admit a KSS vector of infinite kind. In all these 
papers the equation of state has not been specified. 

In recent papers, Maeda et al. [10,11] investigated the KSS vector of 
the second kind in the tilted case. They assumed the perfect fluid spacetime 
obeying a relativistic polytropic equation of state. Further, they assumed two 
kinds of polytropic equation of state in GR and showed that such spacetimes 
must be vacuum in both cases. They studied the case in which a KSS vector 
is not only tilted to the fluid flow but also parallel or orthogonal. 

Daud and Ziad [12] have found the homotheties of spherically symetric 
spacetimes admitting maximal isometry groups larger than SO (3) along with 
their metrics. They have used the homothety equations without imposing 
any restriction on the stress-energy tensor. It would be worth interesting to 
explore the physical properties of these solutions. In this paper, we would 
evaluate the kinematical properties of these self-similar solutions of the first 
kind. We shall discuss the physical properties, such as acceleration, rotation, 
expansion, shear, shear invariant and expansion rate both in co-moving and 
non co-moving (only in radial direction) coordinates. Further, we would look 
for the singularities of these solutions. 

The paper has been organised as follows. In the next section, we shall 
write all such solutions. In section 3, we shall discuss the physical properties 
of these solutions of the first kind both in co-moving and non co-moving 
coordinates. In section 4, singularities of these solutions will be explored. 
Finally, in the last section, we shall conclude the results. 



2 Self-Similar Spherically Symmetric Solutions 

The line element of general spherically symmetric spacetime is given by [13] 

ds 2 = e v{t ' r) dt 2 - e x ^dr 2 - e^dtt 2 , (4) 

where dVt 2 = dQ 2 + sin 2 6d(j) 2 . 

Daud and Ziad [12] have solved homothetic equations and found the dif- 
ferent self-similar solutions of the first kind. There are two classes of such 
solutions, one of which admit 5 HVs and the second class admit 7 HVs. The 
first class of self-similar solutions which admit 5 HVs is given by the following 
three different metrics 

where a is an arbitrary constant but cannot be 1. 

ds 2 = dt 2 -t 2 ^dr 2 -t 2 dn 2 , (6) 

ds 2 = dt 2 -dr 2 -(3 2 {t + ar) 2 dQ 2 , (7) 

where (3 is an arbitrary constant. The second class of self-similar solutions 
which admit 7 HVs is given by the line element 

ds 2 = dt 2 - e x{t) [dr 2 + £ 2 (m,r)dfi 2 ], (8) 

where E = asinh(r/a), r and sin(r/o) for m = 1/a 2 , 0, — 1/a 2 (a being an 
arbitrary constant) respectively subject to the following constraint 

2me~ x + x^0, 

(dot denotes differentiation with respect to t) and 

e x = 2(to + 7 ) 2 , a 2 ^^, 

Or 

= 5(£- 7 r, a^O, 

= (to + 7 ) 2 , «V--, (9) 

a z 

where 7 is an arbitrary constant. Another self-similar solution which also 
admits 7 HVs is given by the following metric 

ds 2 = e y(r) [dt 2 - a 2 cosh(-)dn 2 ] - dr 2 (10) 



subject to the constraint 

2e~ y + a 2 y" ^ 0, 

(prime indicates derivative with respect to r) and 

e y = (ra + [3) 2 , «Vt (11) 

3 Kinematics of the Velocity Field 

In this section, we shall discuss some of the kinematical properties of the 
self-similar solutions given by Eqs.(5)-(8) and (10) both in co-moving and 
non co-moving coordinates. The kinematical properties [13] can be listed as 
follows. The acceleration is defined by 

Ua = U a - b U b . (12) 

The rotation is given by 

Uab = U[ a . b ] +U[ a U b ]. (13) 

The expansion scalar, which determines the volume behaviour of the fluid, 
is defined by 

= u%. (14) 

The shear tensor, which provides the distortion arising in fluid flow leaving 
the volume invariant, can be found by 

1 

Cab = U( a -,b) + U( a U b ) ~ ~Qh ab , (15) 

The shear invariant is given by 

a = a ab a ab . (16) 

The rate of change of expansion with respect to proper time is given by 
Raychaudhuri's equation [14] 

^ = -\Q 2 - a ab u ah + uo ab u a u b - R ab u a u b . (17) 

ar 3 

Now we discuss the properties of these solutions both in co-moving and also 
in non co-moving coordinates. 



3.1 Kinematic Properties in Co-Moving Coordinates 

First we evaluate the kinematical properties of the self-similar solutions in 
the co-moving coordinates. 

For the first solution given by Eq.(5), the acceleration and expansion are 
zero but one of the rotation components exists given by 

I — a , 

u 01 = —. (18) 



The only non-zero component of the shear is 



001 — —Uq\ 


and consequently the shear invariant becomes 


a- {l ~ af 





(19) 



(20) 

By using Raychaudhuri equation, the rate of change of expansion turns out 

to be 

dQ a — I . . 

*=— ■ < 21 » 

For the second self-similar solution given by Eq.(6), the acceleration, rotation 
and expansion are zero. The non-zero shear components are 

on = 2(a - l)^ 1 - 2 *) , a 22 = -2t, a 33 = a 22 sin 2 6. (22) 

The shear invariant becomes 

a=l[(a-l) 2 + 2]. (23) 

The expansion rate is given by Raychaudhuri equation as follows 

^ = --W-9a + 12). (24) 

Now we evaluate the above quantities for the self-similar solution given by 
Eq.(7). It is easy to see that the acceleration, rotation and expansion turn 
out to be zero. The non-zero components of the shear are 

a 22 = -2{3 2 {t + ar) , a 33 = a 22 sin 2 9. (25) 



The shear invariant takes the following form 



o = 7 ^. (26) 

(t + ar) 2 y ' 

The rate of change of expansion is 

* = -■ < 27) 

Now we calculate the above kinematical quantities for another self-similar 
solutions given by Eq.(8). For this class, the acceleration, rotation and ex- 
pansion are zero. The non-zero components of the shear take the form 

On = -x t e x , 

022 = -x t e x £ 2 (ra,r), 

033 = o- 2 2sm 2 9. (28) 

The shear invariant becomes 

a = 3x t 2 . (29) 

The expansion rate turns out to be 

d ^- = -(2x tt -?>x t 2 ). (30) 

Finally, we explore the above quantities for the self-similar solution given by 
Eq.(10). For this solution, the expansion is zero while the acceleration and 
rotation are given by respectively 

ii = ~% (31) 

^oi = 2/re2. (32) 

The non-zero components of the shear are 

001 — _ ^01 5 

°"22 = — 2acosh(-) sinh(-)e2, 
a a 

o"33 = o"22sin 2 #. (33) 



The shear invariant becomes 

a = -y r 2 + 8 ^7^ . (34) 

The rate of change of expansion is given by 



tanh(-' 
a 2 ey 



dB 2a 2 y„e y + 3a 2 y r V - 8 , . 

— = -<r + — — — r^r^ 1 • 35 

dr 4a 2 ev v ' 

3.2 Kinematic Properties in Non Co-Moving Coordi- 
nates 

Here we discuss the kinematical properties of the self-similar solutions in the 
non co-moving coordinates only in radial direction. 

For the first solution, the acceleration components turn out to be 

m = — — , Ml = . (36) 

The rotation component is given by 

u; i = Uq (37) 

and the expansion becomes 

= ^. (38) 

r 

The non-zero components of the shear are 

a 00 = 4(1 - a)r^- 2a \ em = ~, 

Sr 

(a-9)r 2/1 2(4a-3) , , 

0"22 = ~ , (733 = (722 Sin 0, Ooi = — (39) 



and the shear invariant is 



2(45 + 13a 2 - 38a) 
~3? 



a = —2 "■ ( 4 °) 



The expansion rate turns out to be 



d& 1 . 9 . a — 1 , . 

-d7 = ~3? {27a ~ 76a + 93) + ~T- (41) 



For the second solution, the acceleration components become 

a — 1 . (I — a) . . 

Uo = ~T' Ul = ~l^- (42) 

The rotation component is given by 

^01 = -ill (43) 

and the expansion becomes 

e = fc^. (44) 

The non-zero components of the shear are 

a 00 = 2ii , <r 11 = l(5a-3)t<- 1 - 2a \ 

2a 
a 22 = -2t, a 33 = a 22 sm 2 9, cr 01 = - — . (45) 



The shear invariant becomes 



4 ,25 2 16 



a = ^Y a i a + 4) - (46) 

The rate of change of expansion will become 

£» = -J-( 112a . - 219a + 171) + (Ll^ - (1 - ")(" - ") . (47) 

For the self-similar solution given by Eq.(7), the acceleration, rotation and 
expansion are zero. The non-zero components of the shear are 

a 22 = 2[3 2 (t + or) (a - 1), cr 33 = a 22 sin 2 6. (48) 

The shear invariant becomes 

8(1 -a) 2 . . 

a = -± £=-■ 49 

(t + ar) 2 v J 

The expansion rate can be found by using Raychaudhuri equation as given 

below 

dQ . . 

17 = ""■ < 50 » 
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Now we consider the metric given by Eq.(8). In this case, the acceleration 
components turn out to be 

Xt . Xt x . . 

M o = -y, u x = —ei. (51) 

The only non-zero rotation component is given by 

^01 = -ill (52) 

and the expansion becomes 

= |i-2e-f|. (53) 

The non-zero components of the shear are 

o"oo = -x t , (Tn = -2(e s + ^r e ^), °"22 = --e*£E r , 

0-33 = o"22 sin 2 0, (7oi = x t e2 + — ^ (54) 

and consequently the shear invariant becomes 

The rate of change of expansion will be 

dr 2 v 2 ; 2 v ; 3S v ' 

2e 2__ -| 4. (56) 



9e*E 2 S 

Finally, we evaluate the kinematical properties of the self-similar solution 
given Eq.(lO). For this solution, the acceleration components become 

u = —ei, Ul = - — . (57) 

The non-zero rotation component will be 

woi = «o (58) 

10 



and the 


expansion is 


given by 










= A 

aez 


tanh(— ) — 
a 


3 

- 2 Vr 


The non 


-zero shear components are 






o"oo = 


2y r e y , o-ii = 


4tanh(|) 

oaez 


i 




0"22 = 


— acosh(— )[e 


y t 

2(2sinh(-) 


+ y 



(59) 



a a a 

+ -cosh(-)(— s -tanh(-)--j/ r )], 

6 a ae? a z 

o"33 = o"22 sin 2 0, 

y 1 t 

o"oi = -y r e 5 - — tanh(-). (60) 

KjCl Ob 

The shear invariant is 

%2 — -^ l / " \ i 2 / \ 

r y r e 2tanh(-) — tanh (-) 

3a V 3a 2 ey V 

2e~ y t y t 

+ ~ r27Tv(2sinh(-)+|/ r e2acosh(-)) 2 

2e- 2 \ 2 , .t, 3 N2 4e"^ ,2 , ,t, 

+ — 7T- (— tanht-)--^) + - — _(_ s -tanh(-) 

9 ae2 a 2 3a cosh(-) ae2 a 

3 £ t 

— -?/ r )(2sinh(— ) + y r e2acosh(-)). (61) 

The expansion rate is given by 

dQ 1 2 £ 3y r 2 3 3 2 

x tanh(-) — ) + -y rr + -y r 



rfr 3 aet V 2 y 2" 4* 

(2y„.aV + 3y r 2 e y a 2 - 8) - ^ - a. (62) 



4a 2 e^ 

4 Singularities 

In this section, we shall discuss the singularities of the self-similar solutions. 
The Kretschmann scalar is defined by 

K = R abcd R abcd , (63) 

11 



where R a bcd is the Riemann tensor. For the first solution given by Eq.(5), it 
reduces to 

K = ^(l -a)[(l-a)(a 2 + 1) + 1]. (64) 

It is obvious that K diverges at the point r = 0. Thus the solution is singular 
at r = 0. 

For the second solution given by Eq.(6), the Kretschmann scalar turns 
out to be 

K = i[(a 2 + 2)(l-a) 2 + 4]. (65) 

and hence the metric is scalar polynomial singular along t = 0. 
The Kretschmann scalar for the third solution takes the form 

K = — — -\2t 2 + 2atr + a 2 r 2 + 2a(3 2 tr + a 2 (3 2 r 2 

p 4 (t + ra) 8 

- a 2 (3 2 t 2 -2a 3 (3 2 tr-a 4 (3 2 r 2 ]. (66) 

We see that K diverges only when t — 0, r = and consequently this solution 
is singular at these points simultaneously. 

For the class of solutions given by Eq.(9), K becomes 

K = 2[\(x tt + \x?f + -^(E„. - \x t \'Hf 

+ ^(1 + ^W-SVn (67) 

This is singular only when X = 0, i.e., when r = 0. Finally, for the last 
solution, the Kretschmann scalar reduces to 

K = ^Vrr + \yr 2 ? + ^-\yr 2 a 2 ey) 2 



+ 



1 + sinh 2 (t/a) - \y 2 a 2 e v cosh 2 (t/a)) 2 ]. (68) 



a 4 e 2 y C osh 4 (t/a) v 4 



This does not provide any singularity and hence the last solution is singularity 
free. 



5 Conclusion 

In GR, the term self-similarity can be used in two different ways. One is 
for the properties of spacetimes and the other is for the properties of matter 

12 



fields. In general, these are not equivalent. Although homothetic solutions 
can contain many interesting matter fields, but those compatible with homo- 
thety are not so common. In this paper, we have considered some self-similar 
solutions of the first kind and have evaluated some physical properties of these 
solutions both in co-moving and also in non co-moving coordinates. We also 
have checked singularities of each solution (if any). 

First we discuss the physical properties in co-moving coordinates. For 
the first self-similar solution, the acceleration and the expansion turn out to 
be zero while the remaining physical quantities exist and are finite. It can be 
seen from Eqs.(18)-(20) that these quantities cannot be made zero as a ^ 1. 
Further all these quantities become infinite at r = and also the expansion 
rate can be positive or negative according to the choice of a. 

For the second solution, the acceleration, rotation and expansion are zero 
while the rest of the quantities do not vanish as a ^ 1. We see from Eq.(24) 
that the expansion rate does provide positive/negative value according to the 
choice of a. Also, we see that shear invariant and expansion rate turn out to 
be infinite for t = 0. The third solution also gives acceleration, rotation and 
expansion to be zero. In this case the expansion rate is always negative and 
become infinite only at t — 0, r = 0. For the fourth solution again accel- 
eration, rotation and expansion turn out to be zero. Further, the expansion 
rate is zero only if x = constant otherwise it is finite. For the last self-similar 
solution, only expansion is zero while the remaining quantities are non-zero 
and finite. 

In non co-moving coordinates, the kinematical quantities turn out to be 
non- zero for all the solutions except for the solution given by Eq.(7). The 
reason is that we are also considering velocity in radial direction in addition to 
the temporal component. Also, we see that the results are more complicated 
in the non co-moving coordinates as compared to the co-moving coordinates. 

Finally, we discuss the singularities of these solutions. We see that the 
first solution becomes singular at r = while the second solution is singular 
at t — 0. The third solution turns out to be singular only if t — 0, r = 
whereas the fourth solution is singular at r = 0. The last solution turns out 
to be singularity free solution. 
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